PROCEEDINGS 

°F SCIENCE 



o 

< 

(N 



Discovery of VHE gamma-rays from the vicinity of 
tlie shell-type SNR G31 8.2+0.1 with H.E.S.S. 



< 
O 

6 



P. Hofverberg", R.C.G. Chaves", A. Fiasson^, K. Kosack'^, J. Mehaulf^ and E. de Ona 
Wilhelmi ' for the H.E.S.S. Collaboration 

" MPIK Heidelberg, Germany 

^ CNRS/IN2P3, Universite de Savoie, Annecy-le-Vieux, France 
' IRFU/DSM/CEA, CE Saclay, France 
CNRS/IN2P3, Universite de Montpellier 2, Montpellier, France 
E-mail: [petter . hof verbergBmpi-hd . mpg . de 



> 



in 
o 



The on-going H.E.S.S. Galactic Plane Survey continues to reveal new sources of VHE 7-rays. In 
particular, recent re-observations of the region around the shell- type supernova remnant (SNR) 
G318.2H-0.1 have resulted in the discovery of statistically-significant very-high-energy (VHE) 7- 
ray emission from an extended region. Although the source remains unidentified, archival obser- 
vations of '^CO in the region provide an opportunity to investigate a potential SNR/molecular 
cloud interaction. The morphological properties of this newly-discovered VHE 7-ray source 
HESS J1457— 593 are presented and discussed in light of the multi-wavelength data available. 
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1. Introduction 

Supernova remnants (SNR) are thought to be responsible for the acceleration of cosmic rays 
(CRs) up to energies around the "knee" (~10^^ eV) of the CR spectrum, a statement which is 
backed up by both theoretical arguments and experimental evidence. Theoretically, it is well es- 
tablished that supernova explosions release just the right amount of energy into the interstellar 
medium to account for the energy budget of the CRs, assuming they can convert on the order of 
^■10% of their energy into kinetic energy of CRs Furthennore, models exist that can explain 
how CRs can be accelerated through diffusive shock acceleration at the shock front up to energies 
approaching the knee. Experimental evidence comes from primarily X-ray and very-high-energy 
(VHE; E > 0.1 TeV) y-ray studies of young shell-type SNRs such as SN 1006 [@][|], which have 
found evidence for non-thermal populations of CRs extending up to ~ 100 TeV. 

y-rays can be produced in both hadronic interactions (with subsequent pion decay) and from 
energetic electrons that inverse Compton scatter off background photon fields. An ambiguity there- 
fore exists in the responsible radiating particle population when detecting y-rays from SNRs, and 
determining the hadronic or leptonic nature of the accelerating mechanism remains a key issue. 
Since neutrino detectors have not yet proven to be sensitive enough to detect neutrinos from as- 
trophysical sources, no direct indicator of either process currently exists. However, disentangling 
the radiation processes can be facilitated by studying SNR y-ray spectra extending well beyond 
10 TeV, where inverse-Compton spectra tend to cut off due to strong radiative energy losses of the 
parent electron populations, or by studying SNRs interacting with dense molecular clouds where 
the target material in the cloud could boost the y-ray component from a hadronic mechanism. 

For the reasons outlined above, the SNR G3 18.2-1-0.1 may be a prime target for studying VHE 
y-ray emission and the radiation mechanism that generates the emission. G3 18.2-1-0.1 is a large 
diameter (40//x35//) shell-type SNR discovered in the Molonglo Observatory Synthesis Telescope 
(MOST) survey [Q] and is characterized by two non-thermal filaments in the Northwest and South- 
east which form two sections of a shell, and a central region of thermal emission which corresponds 
to a HIT region. 

In this paper, the discovery of the VHE source HESS J1457— 593 is presented and discussed in 
light of multiwavelength data that suggest that the emission originates from an interaction between 
SNR G3 18.2-1-0. 1 and a molecular cloud. 

2. H.E.S.S. observations and Analysis 

2.1 The H.E.S.S. Telescope Array 

The High Energy Stereoscopic System (H.E.S.S.) is an array of Cherenkov telescopes for VHE 
y-ray astronomy by observing Cherenkov light emitted from y-ray induced extensive air showers 
(EASs). H.E.S.S. is comprised of four identical 12 m diameter telescopes arranged in a square 
with 120 m sides and is located in the Khomas highlands of Namibia in the southern hemisphere 
(23°16'17"S) at a height of 1800 m above sea level. The array has been fully operational since 
2004 [Q]. Each telescope is equipped with an alt-az mount and a reflector of Davis-Cotton design 
with a total mirror area of 107 m^ and a focal length of 15 m [^. A photomultiplier camera with 
960 pixels (each of 0.16° size) provides a total field-of-view of 5° which makes H.E.S.S. well 
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suited for observing extended sources and for performing survey observations. The H.E.S.S. array 
uses a stereoscopic trigger which drastically reduces the background from single muons and the 
night sky, keeping dead-time low and allowing a low energy threshold. In addition, employing the 
stereoscopic detection technique, an angular resolution of ~0. 1° and an energy resolution of ~15% 
is achieved. H.E.S.S. 's unprecedented sensitivity to 7-rays permits it to detect a point source with 
a flux of 1% of the Crab Nebula at a significance of 5 a in only ~25 h of observations [||]. 

2.2 Data and Analysis Methods 

The region of interest was observed with the H.E.S.S. telescope array between 2004 and 2010. 
The data-set is primarily comprised not of dedicated observations of the VHE 7-ray source pre- 
sented in this paper but rather of survey observations of the region and from observations of nearby 
sources. After standard quality selection to remove observations taken during unfavorable 
weather conditions or affected by hardware-related problems, the live-time amounts to ~102 h. 
The effective live-time at the source in question is however significantly lower due to the large 
average offset of the pointings relative to the source (1.6°). The data was taken in a series of runs 
of typical duration '^28 min at zenith angles between 35° to 55°, with a mean of 38°. 

The data-set was analyzed using the Hillas second moment method ^ employed in the H.E.S.S. 
standard analysis [Q]. To separate 7-ray event candidates from CR-like events, hard cuts were used, 
where a minimum of 200 photoelectrons is required in each recorded EAS image. Compared to 
standard cuts, which requires a minimum of 80 photoelectrons, this results in a narrower PSF and 
an improved background rejection, but also in an increased energy threshold. The time-dependent 
optical response of the telescopes, due to mirror deterioration, was calibrated using the known 
amount of Cherenkov light from single muons passing close to the telescope 

For the generation of 2D images, the ring background method [||] was used, where the back- 
ground at each point in the sky was calculated within an annulus centered on the point in question, 
with an inner radius of 0.7° and a area factor between the OFF and ON region of 7. The sta- 
tistical significance of the images were derived from the number of off-source (background) and 
on-source events, following the likelihood ratio procedure in [^. All results presented in this pa- 
per were cross-checked with the Model2D analysis [^] which is based on a model of Cherenkov 
image parameters and also utilizes an independent calibration of the raw data. Consistent results 
were obtained with the two analysis techniques. 

3. Results 

Figure |] shows the gamma-ray excess counts in a 2° x 2° region around the source of interest. 
The image has been smoothed with a Gaussian kernel with standard deviation 0.12°. The detection 
significance of the source amounts to 8.9 a using an integration radius of 0.4°. Table 1 gives a 
summary of the event statistics from the analysis. The source centroid and extension were deter- 
mined by fitting an asymmetric 2D Gaussian, convolved with the H.E.S.S. point-spread function, 
to the unsmoothed 7-ray excess map. The best-fit centroid is at aj2ooo = 14'^57™46'* it 33** and 

5/2000 = -59°28' ±5' (1 ~ 318.3°, b 0.4°), giving a ;t;Vndf = 778.7/778, and the source is 

thus named HESS J1457— 593. The best-fit position is indicated in Figure |l| by a white cross, the 
extension of which corresponds to the statistical error of the fit (quoted above). The best-fit source 
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Galactic Longitude (deg) 

Figure 1: Image of the VHE 7-ray excess (in units of 7-rays arcmin^^) of HESS J1457— 593, corrected for 
exposure and smoothed with a 2D Gaussian with a width of a = 0.12°. The color scale is set so that the 
blue/red color transition occurs at the level of statistically significant features (5(7). The black and white 
cross indicates the best-fit position of the source centroid, and its size represent the statistical errors of the 
fit. The fitted extension of the source is indicated by the white dashed ellipse. The inset in the lower right 
corner shows the point spread function for this particular data-set, smoothed in an identical fashion as the 
main figure. The green circle shows the approximate extent of the SNR G318.2H-0.1 in radio. 

extension is amajor = 0.31° ±0.07°, aminor = 0.17° ±0.05, at an angle = 67° ± 15°. It should be 
noted that due to the non-Gaussian morphology of the source, this procedure only gives a rough 
estimate of the source centroid and extension. The north part of HESS J1457— 593 overlaps the 
southern part of the SNR G3 18.2-1-0.1, whose approximate extension in the radio is indicated in 
Figure 1 by a green circle. 

Source Name o;j2ooo 5j2oo() Non Noff a Excess Significance [a] 

HESSJ1457-593 14''57"'46'' -59°28' 5034 22547 019 659 8^9 

Table 1: HESS J1457— 593 source centroid (columns 2 and 3), event statistics (columns 4-7) and peak 
significance (column 8). 



4 



Discovery ofVHE gamma-rays from the vicinity of the shell-type SNR G318. 2+0.1 with H.E.S.S. 



4. Observations of molecular clouds 



In searching for molecular cloud counterparts, ^^CO data from the Dame survey [11] was 
analyzed. A square box, 1° x 1° large, centered on the best fit position of HESS J1457— 593 was 
used as an integration region to search for evidence of molecular clouds in the region of the H.E.S.S. 
source. Integrating ^^CO over the full Local Standard of Rest velocity (Vlsr) range, a strong peak 
was found in the velocity interval —54.5 km/s to —32.5 km/s centered at "Vlsr = —42.0 km/s. A 
map of the ^^CO line emission integrated in this velocity interval is shown in Figure The color 
transition is chosen so that the blue-red transition start at the boundaries of the molecular cloud, 
derived using the definition of molecular cloud boundaries suggested by [|l2|]. The cloud complex 
coincident with the H.E.S.S. source measures roughly 1.8° x 1.1° and fully contains the VHE y-ray 
emission (indicated by black contours). 

The distance to the molecular cloud is estimated assuming the Galactic rotation curve from 
[13] with a rotation speed of Vq = 217 km/s and a Galactocentric distance of the sun of Rq = 



8.5 kpc. A velocity distance of —42.0 km/s then corresponds to a near/far kinematic distance of 
(3.5 ±0.2) kpc and (9.2 ±0.2) kpc respectively, where the errors are estimated by comparing the 



results above with the results obtained using a Galactic rotation model from [14] . 

Assuming a relationship between the velocity-integrated CO intensity Wco and the molecular 
hydrogen column density N(H2) as N(H2)/Wco = 1.8 x lO^*' [11], the cloud properties given in 



Table 2 are derived for the near kinematic distance solution. Using the far distance solution, the 
molecular cloud mass would be in the order of 10^ M©. Such massive clouds have not been 
observed outside the Galactic centre and the near distance solution is therefore adopted throughout 
this paper. In the above calculations, it is assumed that the cloud extension in the line-of-sight 
is equal to the average cloud extension perpendicular to the line-of-sight. The molecular cloud 
coincident with the H.E.S.S. source has the characteristics of a typical Giant Molecular Cloud 



(GMC) with a mass on the order of 10 Mg, and an extension of tens of parsec [15]. 



Position Angular Extension Average Extension Mass Density 

(l,b)=~(318.4°,-0.5°) 1.8° X 1.1° ~80 pc ~3 x lO^Mg ~40 cm^^ 

Table 2: Derived quantities of the GMC coincident with the source HESS J1457— 593. The first column 
gives the geometrical centre position of the cloud, the second its angular extension (in Galactic longitude, 
latitude). Columns three to five gives the average (physical) extension, the mass and the density of the cloud, 
all assuming a distance to the cloud of 3.5 kpc. 



5. Discussion and Conclusions 

An extended source of VHE y-ray emission, HESS J1457-593, has been found in the vicinity 
of the SNR G3 18.2-1-0.1. As shown in Figure |, the TeV y-ray emission is partially overlapping 
the non-thermal, southern rim of the SNR and has an elongated morphology. A GMC complex is 
positionally coincident with the TeV y-ray emission and overlapping the southern rim of the SNR, 

'The distance uncertainty arising from the velocity resolution of the '^CO survey, 1.3 km/s, is negligible. 
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Figure 2: Map of the '^CO (1:1-^-0) line emission around HESSJ1457-593 integrated between -54.5 km/s 
and —32.5 km/s. Overlaid in turquoise are the excess count contours of HESS J1457— 593 (corresponding 
to 0.40, 0.50 and 0.60 7-rays arcmin"^). In green are the radio contours (843 MHz) of the SNR G318.2H-0.1. 
The contours correspond to intensities of 2, 6 and 10 mJy /beam. The colorbar units are in Kkm/s. 



suggesting a region of cosmic -ray enhancement and providing an excellent opportunity to inves- 
tigate a putative SNR-GMC interaction and to study cosmic ray propagation through a molecular 
cloud. 

Assuming the GMC and the SNR G318.2-H0.1 are spatially associated, the distance to the 
SNR, unknown until now, is given. This is a key parameter of SNRs and allows an estimation of 
its physical size and age. At a distance of 3.5 kpc and with an angular size of (40//x35//), the SNR 
would have a physical diameter of roughly 40 pc. A Sedov-Taylor model of SNR evolution [16] is 
then used to estimate its age directly from the SNR diameter and the ambient density which can be 
assumed to be 1 cm^^, the average for Galactic SNRs. This gives an age of SNR G318. 2-1-0.1 of 
roughly 8000 years^. SNR G3 18. 2-1-0.1 would then be in the late Sedov phase and could thus be an 
efficient particle accelerator. 
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